Quantum algorithms without initializing the auxiliary qubits 
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In this Letter, we construct the quantum algorithms for the Simon problem and the period-finding 
problem, which do not require initializing the auxiliary qubits involved in the process of functional 
evaluation but are as efficient as the original algorithms. In these quantum algorithms, one can use 
any arbitrarily mixed state as the auxiliary qubits, and furthermore can recover the state of the 
auxiliary qubits to the original one after completing the computations. Since the recovered state 
can be employed in any other computations, we obtain that a single preparation of the auxiliary 
qubits in an arbitrarily mixed state is sufficient to implement the iterative procedure in the Simon 
algorithm or the period-finding algorithm. 
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Quantum computational algorithms can be executed 
in parallel on superpositions of exponentially many input 
states, and their outcomes can be properly measured by 
virtue of quantum interference. These enable exponential 
speedups in the solutions of certain problems, and allow 
one to distinguish between the quantum computational 
complexity classes and the classical ones pj 0,13, [3- 13- E3- 

HHSE3 

It is assumed that most quantum algorithms require 
some initialization at start-up, which is to prepare a cer- 
tain pure state as an initial state. However, in experi- 
mentally realizable proposals for the implementation of 
quantum algorithms, it may be technically difficult to 
prepare the initial pure state. Especially, the nuclear 
magnetic resonance (NMR) system is typically applied 
to physical systems in equilibrium at room temperature. 
This means that the initial state of the spins is nearly 
completely random, that is, it is difficult to prepare pure 
quantum states of nuclear spins in NMR systems. Hence, 
it would be interesting whether quantum algorithms can 
be efficiently performed even though the initial states or 
some parts of them are not in a specific pure state. If it 
would be possible, then the technical difficulty could be 
settled to a certain extent, and furthermore if the parts of 
the initial state would remain intact even after the com- 
putation, then the parts could be reused in any other 
computations. We call such a quantum algorithm the 
initialization-free quantum algorithm when any quantum 
state can be used as the auxiliary (target) qubits involved 
in the functional evaluation |x) <8> \y) t— > \x) ® \y + f(x)) 
for a given function /, and it can be recovered after the 
computation. 

In the initialization- free quantum algorithms, which 
are implemented without initializing and deforming the 
state of the auxiliary qubits, any qubits (which might 
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contain some other useful information) can be temporar- 
ily used as the auxiliary qubits, and the initial state of 
the auxiliary qubits can be recovered at the end of the 
computation. Thus we can compose the auxiliary qubits 
of any qubits regardless of whether they are entangled 
with others or being used in another computational pro- 
cess. Furthermore, in the case of iterative algorithms, in 
which one needs to perform the algorithm several times 
to solve given problems, the initialization- free quantum 
algorithms can be implemented with the same auxiliary 
qubits repeatedly, while the original iterative algorithms 
require the initial auxiliary qubits of a certain pure state 
at each repetition. 

There have been a few research works related to 
the initialization-free quantum algorithms. Biham et 
al. [Ill IT^ I have generalized Grover's algorithm |(| by al- 
lowing for an arbitrary initial amplitude distribution, and 
have shown that Grover's algorithm (or, a large class of 
Grover-type algorithms) is robust against modest noise in 
the amplitude initialization procedure. Parker and Ple- 
nio |I3| found that one pure qubit and an initial supply of 
log 2 N qubits in an arbitrarily mixed state arc sufficient 
to implement Shor's quantum factoring algorithm [j| effi- 
ciently, where the idea of using one pure qubit and other 
mixed qubits as the initial state was first introduced by 
Knill and Lafiamme [l4|- Their result implies that the 
controlled unitary transformations in Shor's algorithm 
can be implemented without any initialization of the aux- 
iliary qubits, while the auxiliary qubits cannot be left in 
the initial state. Subsequently, Chi et al. 15] have pre- 
sented a quantum algorithm to implement an oracle com- 
puting \x) i — ► e**if(.x)/M \ x ) for / : Z N -> Z M by making 
use of an oracle of the form \x) ® \y) i— > \x) ® \y + f(x)) 
without setting the auxiliary qubits to a definite state be- 
fore the computation, and have shown that generalized 
Deutsch-Jozsa algorithms can be implemented without 
any initialization of the auxiliary qubits when an oracle 
computing Uf : \x) (8 \y) i— > \x) ® \y + f(x)) is employed. 

In this Letter, we deal with two problems, the Si- 
mon problem 4] and the period-finding problem pj, 
which can be solved efficiently by the quantum com- 
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puter, and present the initialization-free quantum algo- 
rithms for these problems. Since most of known "expo- 
nentially fast" applications of the quantum Fourier trans- 
form (QFT) can be considered as a generalization of find- 
ing unknown period of a periodic function, the existence 
of these quantum algorithms implies that any initializa- 
tion of the auxiliary qubits may be unnecessary in many 
quantum algorithms. 

We first recall the Simon problem £| that can be solved 
in polynomial time on a quantum computer but that re- 
quires exponential time on any classical bounded-error 
probabilistic Turing machine if the data is supplied as 
a black box, and we then investigate initialization-free 
techniques for the Simon problem. We note that there 
exists an exact quantum polynomial-time algorithm for 
the Simon problem 9]. However, we here deal with the 
original Simon algorithm which is polynomial-time in the 
expected sense. 

For convenience, we use the following notations. Let 
G = {HQ > ffin) be a group under the binary operation © n , 
which is the bitwise XOR operation. For a subset A of 
G, let \A\ denote the cardinality of A. 

We define a bilinear map G x G — > Z2 by 



x-y = x y © ccij/i © • • • © x„_ij/ n _i 



(1) 



where x = (xq,xi,- ■ ■ ,ar n -i) and y = (2/0,2/1,- ■ ■ ,2/n-i) 
are elements in G for Xj , yj E Z2 , and © is the XOR 
operation, which is the addition modulo 2. This bilinear 
map clearly satisfies the property that (x ©„ y) ■ z = 
x ■ z © y ■ z for x, y and z in G. 
For a subgroup H of G, let 



H 1 - = {x e G : x ■ y = for all y e H) 



(2) 



denote the orthogonal subgroup of H . We remark that 
the quotient group G/H is well-defined since G is an 
abelian group. 

Let / : G — * G be an arbitrary two-to-one map such 
that f(x) — f(y) if and only if x ©„ y € H where H = 
{0, h} is a subgroup of G for some nonzero h £ G. Then 
the Simon problem is to find the subgroup H , that is, to 
determine the value of h. The original Simon algorithm 
is as follows: (i) Prepare |0 n ) © |0 n ). (ii) Apply W„ © 
I, where W n is the n-qubit Walsh-Hadamard transform 
defined as \x) ~ (1/ y/ffi) "E yeG {-l) x - v \y). (Hi) Apply 
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U f . (iv) Apply W„ © T. Then the resulting state is 

i $ ) = i|| E E (-n^i/w). (3) 

' ' yeH 1 - xEG/H 

We measure the first n-qubit state. Then for each 
y 6 H , the probability with which we obtain y as the 
outcome is 

(*|(|l/>(tf|® I)|$) = pa E 1 = T^- ( 4 ) 



xGG/H 



\Gl 



Thus, after expected 0(n) repetitions of this procedure, 
at least n linearly independent values of y can be collected 
so that the nontrivial h* is uniquely determined by solv- 
ing the linear system of equations h* ■ y = 0. Then we 
have h* = h as required. 

Now we present the initialization-free quantum algo- 
rithm for the Simon problem. We first consider the fol- 
lowing quantum algorithm: (i) Prepare an n-qubit state 
in the state |0 n ) as the control qubits and an n-qubit 
state in an arbitrary pure state |^) = X)/c a fel^) as the 
auxiliary qubits. (ii) Apply W n ®I. (Hi) Apply W®, 
where U® : \x) © \y) 1— > \x) © \y ©„ f{x)). (iv) Choose a 
random n-bit string w — (wq,Wi, ■ ■ ■ ,w n -i) and apply 
S w = er™° © af 1 © • • • © cr™" _1 on the n-qubit auxiliary 
qubits, that is, apply X © S w . (v) Apply Uf. (vi) Apply 
T © S w . (vii) Apply W n © T. Then the resulting state 
becomes 



w\ E E (-ir'<->(-ir 



I*)- (5) 



We now measure the first n-qubit state. Then for each 
y £ H- 1 , the probability with which we obtain y as the 
measurement outcome is 



P w {y) 



|GP 



E (-ir /(x) (-i; 

x<£G/H 



x-y 



(6) 



Hence, the expected probability of obtaining y for ran- 
domly chosen w is 



igi E p »(y) - | G | 3 E 



ibEG 



x<£G/H 



,4 E (E(-ir (/( ^ 

1 1 x,x'£G/H \w£G 



)ffin/(x')) j r_l)(^ena: / ) , I/ < 



(7) 



Since / is one-to-one on G/H, that is, /(x) 7^ /(^') if and only if a; 7^ a;', the inner summation in J7J always 
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vanishes for f(x) ® n f(x') ^ 0, and the summation is 
\G\ for f(x) ©„ f(x') = 0. Thus, for each y € H ± , the 
expected probability (JZJ) becomes 2/|G|. 

For any auxiliary qubits of the state ps = 
£fcPfcl**><*k|» we let p = |0")(0™| ® pb- Then the su- 
peroperator A, which maps p to the quantum state 



1 



(8) 
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where A,,, = (W„ ®T)(I® S w )Uf{X <g> S w )Uf{W n ® X), 
performs the initialization-free Simon algorithm, since if 
the first n-qubit state in A(p) is measured, then it follows 
from J7J that the probability to obtain y as the measure- 
ment outcome is 

tr[{\y){y\®l)k{p)] = ± £ P w (y) = A, ( 9 ) 



tuEG 



|G| 



which is the same probability as that of the original 
Simon algorithm. Furthermore, when y is obtained 
as the measurement outcome, the resulting state after 
the measurement becomes <8 Pb- Therefore, this 

initialization-free quantum algorithm can efficiently solve 
the Simon problem. 

Similarly, we can present the initialization-free quan- 
tum algorithm for the period-finding problem. We first 
review the original quantum algorithm for the period- 
finding problem, and then present the initialization-free 
period-finding algorithm, which can be considered as a 
generalization of the initialization-free Simon algorithm. 

Let / : Z2" — > Z2™ be a periodic function with 
an unknown period T, that is, f(x) — f(x + kT) for 
< k < [2 n /T\ (or, < k < [2 n /T\ + 1). Then 
the period-finding problem is to find T. Classically, this 
problem is known to be hard in the sense that no clas- 
sical algorithms which can find T in polynomial time 
have been found. However, there exists a polynomial- 
time quantum algorithm for the period finding 5] , which 
is as follows: Let N = 2™. (i) Prepare |0") <g> j0 m ). (ii) 
Apply T®X, where T is the Af-dimensional QFT defined 

as \x) 1 * (l/VN)Y,y=oZ 2nixv/N \y)- R Apply W/. (iv) 
Apply T ® X. Then the resulting state becomes 

iEEE <?^^\y) 9 \f(x)) (10) 
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y = x = j=0 



where A x = \_N/T\ or \N/T\ + 1. Now we measure the 
first n-qubit state, and then the probability of obtaining 
y as a measurement outcome is 



1 

p (y) = E 



E< 



,2myjT/N 



(11) 



We note that there are precisely T values of y in 
{0, ,N-1} satisfying 



- < yT (mod N) < -, 



(12) 



and for each y satisfying l|12|) , the probability of obtaining 
such y can be bounded asymptotically, 



(13) 



Thus, with probability at least 4/7r 2 , the measured value 
of y satisfies the inequalities (I12f> . that is, y satisfies the 
following inequalities: 



k 
T 



1 y k 
< — < — 

2N ~ N ~ T 



1 

27V' 



or equivalently 



JJ_ 

N 



< 



2N 



(14) 



(15) 



with k randomly chosen in {0, 1, • • • , T — 1} depending 
on the measurement outcome. Therefore, for sufficiently 
small T with respect to N, the value k/T can be effi- 
ciently extracted from the measured y/N by the contin- 
ued fraction method. Since k and T may be relatively 
prime with high probability, we can get the period T in 
polynomial time with respect to log N . 

The initialization-free quantum algorithm for the 
period-finding problem can be presented by the proce- 
dure similar to the initialization-fee Simon algorithm. In- 
stead of S w = a™ ® af 1 ® • • • ® cr™" -1 for a randomly 
chosen n-bit string w in JSJ), we employ an m-qubit uni- 
tary operation 1l w for a randomly chosen m-bit string w 
defined as \y) i-> e 27riwy ^ M \-y), where M = 2 m_ We pr0 _ 
ceed with the following quantum algorithm: (i) Prepare 
an n-qubit state in the state |0 ra ) and an m-qubit state in 
an arbitrary pure state |^) = £ fc at\k). (ii) Apply T®T. 
(iii) Apply Uf. (iv) Choose a random m-bit string w and 
apply IA W on the m-qubit state of the auxiliary qubits, 
that is, apply I&U W . (v) Apply U / . (vi) Apply T®U W . 
(vii) Apply T <E> 1. Then the resulting state becomes 



N-l /T-11,-1 



Jf E E E e 2 ^ x+jT V N e 2 ™fW/ M \y) ® 



y=0 \a;=0 j=0 



(16) 



Hence, the probability with which we can get \y) as a 
measurement outcome of the first 71-qubit state is 



P w (y) 



1 

iV 2 



T-l A»-l 



2my{x+jT)/N 2iriwf(x)/M 



EE- 



(17) 



By straightforward calculations, we can get the expected 
probability of obtaining y for randomly chosen w, 



^ M-l T-l 

sE p -m = i^E 



M 



tu=0 



z=0 



E^ 

j=0 



2myjT/N 



(18) 
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As in the initialization-free Simon algorithm, for any 
m-qubit state p B = EfcP*|*fc> (*fe|. we let p = |0")(0"|® 
ps, and let the superoperator A be defined as 



M-l 



(19) 



uj=0 



where A,,, = (T ®l){l®U w )U f {l ®U w )U f {T Then 
the superoperator A can perform the period-finding al- 
gorithm efficiently without any initialization on the aux- 
iliary qubits, since the probability of obtaining \y) satis- 
fying JT2J) is 



tx[(\y){y\^l)K{p)} = l^ 



3=0 



2iriyjT/N 



> (20) 



which is the same probability as that of the original 
period- finding algorithm in Furthermore, as in the 

initialization-free Simon algorithm, the resulting state 
after the measurement becomes \y}{y\ <8 Pb when y 
is the measurement outcome. Therefore, there exists 
an initialization-free quantum algorithm which can ef- 
ficiently solve the period-finding problem. 

In conclusion, we have investigated the initialization- 
free quantum algorithms, which do not require any ini- 
tialization of the auxiliary qubits involved in the process 



of functional evaluation, and which recover the initial 
state of the auxiliary qubits after completing the compu- 
tations. We have considered quantum algorithms for the 
Simon problem and the period-finding problem, and have 
presented the initialization-free quantum algorithms for 
the problems, which are as efficient as the original ones. 

The iterative algorithms such as the Simon algorithm 
and the period-finding algorithm demand the storage of 
auxiliary qubits and the extra operations to initialize the 
state of the auxiliary qubits, whenever the procedure re- 
peats. However if one utilizes our initialization-free tech- 
nique then the size of the storage can be reduced and the 
extra operations can be omitted, since the same auxiliary 
qubits can repeatedly be used in our algorithms. Further- 
more, since most known applications of the QFT can be 
considered as a generalization of finding unknown period 
of a periodic function (for example, Shor's factoring algo- 
rithm [f| and Hallgren's more recent algorithm for solv- 
ing Pell's equation |10|). the initialization- free technique 
could be applied to a lot of implementations of quantum 
algorithms. 
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